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Abstract.  This  is  the  continuation  of  the  effort  to  utilize  a  kinetic  moleciilar  model  approach  to  the  problem 
of  turbulence  oriented  computation  of  compressible  flow.  The  Boltzmann  equation  is  employed  in  its  integral 
form  with  the  BGK  model  of  the  collision  term.  This  time,  we  compute  a  periodic  compressible  flow  inside 
a  unit  cube  from  a  random  velocity  field  having  an  energy  spectrum  of  isotropic  type,  uniform  density  and 
temperature.  Results  show  the  development  of  many  small  shock  waves  like  structures  along  with  vortex 
or  eddy  shocklets  scattered  in  the  entire  flow  field.  While  the  energy  spectrum  does  not  change  much  in 
its  pattern  with  time,  as  should  be  the  case,  .some  quantities  like  density  distribution  changes  quickly  to  a 
turbulent  state  from  its  initial  uniform  one.  Some  geometric  properties  of  computed  flow  field  are  derived 
from  its  velocity  deformation  tensor. 


INTRODUCTION 

We  have  been  concerned  with  the  kinetic  molecnlar  model  approach  to  the  problem  of  turbulence  oriented 
computation  of  compressible  flow  and  considered  the  Taylor-Green  type  initial  value  problem  [1],[2].  This  was 
to  see  the  phenomenon  of  the  production  of  small  eddies  from  large  ones,  whose  mechanism  was  considered  to 
be  one  of  the  basic  features  of  turbulent  flow.  Here  we  consider  the  feature  of  developed  turbulent  field.  There 
have  been  some  attempts  to  this  problem  with  use  of  lattice  Boltzmann  equation  [3].  Also  more  recently  an 
interesting  preliminary  study  for  this  problem  with  DSMC  approach  is  reported  [4].  We  compute  a  periodic 
compressible  flow  inside  a  cube  started  from  a  random  velocity  field  with  an  energy  spectrum  having  the  nature 
of  an  isotropic  turbulence.  We  use,  as  in  the  previous  investigations  the  Boltzmann  equation  in  its  integral 
form,  which  is  reduced  for  small  time  step  At  as 

f{c,x  +  cAt,t-\-At)  -  (1) 

where  /  =  f(c,  x,t)  is  the  molecular  velocity  distribution  function  with  the  molecular  velocity  c  and  the  spatial 
coordinate  x,  and  def/dt  denotes  the  rate  of  change  in  molecular  particle  number  owing  to  their  encounters, 
fn  practice,  we  utilize  RGK  model  for  the  term  def/dt,  and  evaluate  the  value  at  x  —  cAt,  which  is  not 
necessarily  a  mesh  point  with  use  of  a  linear  interpolation  from  triangle  (for  20  case)  or  tetrahedron  (for  30 
case)  of  neighboring  mesh  points.  The  consistency  of  the  present  method  with  the  other  approaches  such 
as  with  a  different  collision  model  or  the  Navier-Stokes  CFO  method  has  been  observed  well  in  the  previous 
investigation  [1].  Note  also  both  the  similarity  and  the  difference  between  eq.(l)  and  the  lattice  Boltzmann 
equation  [3]  that  is  a  discrete  model. 


ISOTROPIC  TURBULENT  FIELD 

Consider  here  an  initial  value  problem  to  an  isotropic  turbulent  field  that  is  periodic  to  a  cube  region.  We 
set  for  this  a  random  initial  condition  at  the  time  given  by  a  local  Maxwellian  with  its  flow  velocity  w,,  density 
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FIGURE  1.  Tso-siirface  of  density  (a)  and  vorticity  (b)  distributions  at  t=0.0125  for  A^„=0.05. 


p  and  temperature  T  for  x  =  {x,y,z)  where  the  magnitude  of  u  is  given  by  the  energy  spectnim  E{k)  of 
isotropic  nature  as 

E{h),E  =  2(Ar/ATo)^},  ko  =  10,  k  =  |fe|  for  wave  vector  k,  (2) 

while  its  phase  is  givenin  random,  and  p,  T  are  set  uniform  for  simplicitys  sake.  Here  all  quantities  are 
dimensionless  based  on  the  length  of  an  edge  of  the  cube,  the  initial  uniform  density  and  temperature.  We  use 
Mach  number  M  =4.3  (E  =  0.07)  to  this  initial  field  given  by  the  velocity  to  give  the  maximum  energy  in 
fe-space.  We  have  a  relation  between  the  Reynolds  number  R^.,  the  Knudsen  number  Kn  and  the  Mach  number 
M  as  ReKn  =  SM,  {S  =  8/5i/^  ~  0.9  for  for  a  Maxwell  molecular  model  ). 


RESULTS 


Computational  range  of  »  is  a  cube  region,  n-space  is,  as  usually  been  practiced,  restricted  to  a  finite 
region  |e|  <4  by  setting  outside.  We  set  25x25x25  divisions  for  x,  and  8x8x8  divisions  for  c.  Results  give 
distributions  of  p,  p,  T  and  u  over  c  at  various  time  levels.  These  are  also  utilized  to  derive  some  quantities 
characteristic  to  the  flow  field  such  as  the  energy  spectrnm  E{k),  and  invariant  of  the  deformation  tensor  of 
velocity  field.  Some  representative  examples  from  these  results  are  shown  in  the  following  figures.  Fig.l  shows 
iso-surfaces  of  density  (a)  and  vorticity  (b)  distributions  at  t=0.0125  for  /fn=0.05  so  that  the  initial  Reynolds 
number  i7,e=128  for  M  =4.3.  We  can  see  the  density  distribution  changes  quickly  to  a  turbulent  state  from  its 
initial  uniform  one.  We  see  both  in  (a),  (b)  that  the  development  of  many  small  shock  waves  likes  structures 
combined  with  vortices  and  eddy-  shocklets  [5]  scattered  in  the  entire  flow  field.  Fig. 2  (a)  shows  the  energy 
spectra  at  various  time  levels  for  Rn=0-05.  The  pattern  do  not  change  much  with  time,  as  it  should  be  to  the 
case,  and  they  all  show  the  tendency  of  deviating  of  the  slop  of  curve  from  -5/3  in  contrast  to  the  Taylor-Green 
type  flow  case  [2].  Some  geometric  properties  of  flow  field  [6]  are  studied  by  examine  the  velocity  deformation 
tensor  E  =  {fiij},  where  fijjis  its  components.  The  values  of  its  invariant  P,  Q  and  R  defined  as 

P  =  -det(fij),  Q  =  fill  622  +  fi22fi3.1  +  63.3611  -  6236.32  “  612621  -  61.2fi.21  =  divt?  (3) 

are  computed  for  every  mesh  points  of  x,  at  some  instants  t  and  they  are  plotted  in  planes  (P,Q),  (Q^R)  and 
(P,  R).  Of  these  (Q,  i?,)-plots  at  t=0.0125  are  shown  in  Fig. 2(b).  We  can  notice  there  a  tendency  of  these  points 
concentrating  near  the  line  Q  =  —R,  indicating  the  characteristics  of  the  flow  field  in  compressible  turbulence. 
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(a)  (b) 

FIGURE  2.  (a)  Energy  spectra  and  (b)  spatial  distribution  of  invariant:(Q,  /?)-plot  of  velocity  deformation  tensor  at 
*=0.0125/fN=0.05  for  /<'n=0.05. 
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